Emerging infectious diseases and animal social systems

Emerging infectious diseases threaten a wide diversity of animals, and important questions remain concerning disease emergence in socially structured populations. We developed a spatially explicit simulation model to investigate whether—and under what conditions—disease-related mortality can impact rates of pathogen spread in populations of polygynous groups. Specifically, we investigated whether pathogen-mediated dispersal (PMD) can occur when females disperse after the resident male dies from disease, thus carrying infections to new groups. We also examined the effects of incubation period and virulence, host mortality and rates of background dispersal, and we used the model to investigate the spread of the virus responsible for Ebola hemorrhagic fever, which currently is devastating African ape populations. Output was analyzed using regression trees, which enable exploration of hierarchical and non-linear relationships. Analyses revealed that the incidence of disease in single-male (polygynous) groups was significantly greater for those groups containing an average of more than six females, while the total number of infected hosts in the population was most sensitive to the number of females per group. Thus, as expected, PMD occurs in polygynous groups and its effects increase as harem size (the number of females) increases. Simulation output further indicated that population-level effects of Ebola are likely to differ among multi-male–multi-female chimpanzees and polygynous gorillas, with larger overall numbers of chimpanzees infected, but more gorilla groups becoming infected due to increased dispersal when the resident male dies. Collectively, our results highlight the importance of social system on the spread of disease in wild mammals.     

Simulation of Y-chromosomal haplotype data

The non-recombining nature of the Y-chromosome determines the non-independence of alleles between loci. The evolution of short tandem repeat (STR) loci in the Y-chromosome is the result of different factors such as differential mutation rates, mutation modes, gene conversion, selection and demographic processes. The degree of correlation between loci is dependent on the magnitude of these processes. The simulation of data is a routine tool used for testing hypotheses in population and evolutionary studies. The most basic parameters hitherto used in lineage haplotype simulations are the allele frequency distributions and mutation rates, assuming either full independence or linkage between loci. In this study we introduce use of the Spearman correlation coefficient to estimate the degree of dependence between non-recombining loci. Then, both the interdependence between loci and the allele frequency distributions at multi-allelic loci are incorporated in an algorithm for simulating haplotypes. We illustrate the method using published and unpublished Y-chromosome STR data.     

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W1 and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the O_b atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.     

A Molecular Dynamics Investigation of Chain Conformational Changes In Compressed Bilayers

Abstract

The conformations of the chains constituting the hydrophilic component of alkyl monolayers and bilayers are investigated by performing molecular dynamics atomistic simulations on these systems at different temperatures. Several monitoring techniques are used to reveal the chain conformations, including atom pair radial distribution functions, evolutions of the torsional angles over thousands of timesteps, frequency distributions of the torsionl angles and ‘snapshot’ plots of the atomic configurations. These methods consistently testify to the stability of the trans (fully extended) character of the strain-free alkyl chains up to room temperature. The chains retain much of this conformation even when the layers are compressed by the application of pressure, to which the chains respond by ‘folding’ at the ends attaching them to the substrate planes while maintaining directions which are mainly normal to these planes. A non-zero gap between the layers is also maintained. A pressure of about 50 kbar abruptly causes all motion in the chains to cease, resulting in a highly ordered lattice structure.     

Computational Oral Absorption Simulation for Low‐Solubility Compounds

Bile micelles play an important role in oral absorption of low‐solubility compounds. Bile micelles can affect solubility, dissolution rate, and permeability. For the pH–solubility profile in bile micelles, the Henderson–Hasselbalch equation should be modified to take bile‐micelle partition into account. For the dissolution rate, in the Nernst–Brunner equation, the effective diffusion coefficient in bile‐micelle media should be used instead of the monomer diffusion coefficient. The diffusion coefficient of bile micelles is 8‐ to 18‐fold smaller than that of monomer molecules. For permeability, the effective diffusion coefficient in the unstirred water layer adjacent to the epithelial membrane, and the free fraction at the epithelial membrane surface should be taken into account. The importance of these aspects is demonstrated here using several in vivo and clinical oral‐absorption data of low‐solubility model compounds. Using the theoretical equations, the food effect on oral absorption is further discussed.     

Field studies on the population dynamics of Biomphalaria glabrata, intermediate host of Schistosoma mansoni on the West Indian Island of St. Lucia

Growth curves, calculated for field populations of B. glabrata, were not materially affected by habitat, altitude or season. A mean growth curve was therefore used to estimate the age-frequency distribution of snails in successive field samples. These data permitted the construction of ecological life tables and the estimation of r, the intrinsic rate of natural increase (or decrease) of the different populations. The calculated values of r were inserted in a simple model of unlimited population growth but the resultant curves poorly represented the observed data. A model for unlimited growth was more satisfactory for pond and marsh populations but, apparently, immigration made it less satisfactory for stream and banana drain populations. Nevertheless, r may still be of value in predicting repopulation rates in certain habitats after a mollusciciding which does not kill the entire snail population.     

Trait plasticity alters the range of possible coexistence conditions in a competition–colonisation trade‐off

Most of the classical theory on species coexistence has been based on species‐level competitive trade‐offs. However, it is becoming apparent that plant species display high levels of trait plasticity. The implications of this plasticity are almost completely unknown for most coexistence theory. Here, we model a competition–colonisation trade‐off and incorporate trait plasticity to evaluate its effects on coexistence. Our simulations show that the classic competition–colonisation trade‐off is highly sensitive to environmental circumstances, and coexistence only occurs in narrow ranges of conditions. The inclusion of plasticity, which allows shifts in competitive hierarchies across the landscape, leads to coexistence across a much broader range of competitive and environmental conditions including disturbance levels, the magnitude of competitive differences between species, and landscape spatial patterning. Plasticity also increases the number of species that persist in simulations of multispecies assemblages. Plasticity may generally increase the robustness of coexistence mechanisms and be an important component of scaling coexistence theory to higher diversity communities.     

Application of a photosystem II model for analysis of fluorescence induction curves in the 100 ns to 10 s time domain after excitation with a saturating light pulse

A mathematical model of photosystem II (PSII) events was used to analyze chlorophyll fluorescence transients in the time domain from 100 ns to 10 s after excitation with a saturating 10-ns flash, applied as a part of specialized illumination protocol, using preparations of a thermophilic strain of the unicellular green alga, Chlorella pyrenoidosa Chick (using both intact and diuron-treated cells). Analysis of simulation results has proven that particular attention should be given to flash-induced recombination processes, including nonradiative recombination in PSII, while subsequent charge transfer along the electron transport chain of thylakoid membrane can be adequately described by a single reaction of quinone reoxidation. The PSII model was extended by taking inhibition by diuron of the electron transport in the acceptor side of PSII into account, which allowed simulation of fluorescence induction curves observed in the presence of this inhibitor. The model parameters were determined (stromal pH, rate constants of nonradiative recombination, and the initial reduction state of the quinone pool) which provided adequate simulation of experimentally observed ratios of the maximal and initial fluorescence levels (F _m/F ₀).